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Abstract 



We describe a rearrangement of the standard expansion of the symmetry 
breaking part of the QCD effective Lagrangian that includes into each or- 
der additional terms which in the standard chiral perturbation theory (%PT) 
are relegated to higher orders. The new expansion represents a systematic 
and unambiguous generalization of the standard xPT, and is more likely to 
converge rapidly It provides a consistent framework for a measurement of 
the importance of additional "higher order" terms whose smallness is usually 
assumed but has never been checked. A method of measuring, among other 
quantities, the QCD parameters m{qq) and the quark mass ratio m s /m is 
elaborated in detail. The method is illustrated using various sets of available 
data. Both of these parameters might be considerably smaller than their re- 
spective leading order standard xPT values. The importance of new, more 

accurate, experimental information on low-energy it — tt scattering is stressed. 
PACS numbers: 11.30.Rd,11.40.Fy,13.75.Lb 
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I. INTRODUCTION 



Owing to quark confinement, the connection between QCD correlation functions and 
hadronic observables is far from being straightforward. In the low-energy domain, such a 
connection is described by chiral perturbation theory (%PT) |T|-|3|]. The latter provides a 
complete parametrization (in terms of an effective Lagrangian) of low-energy off-shell cor- 
relation functions of quark bilinears,which should take into account: (i) the normal and 
anomalous Ward identities of chiral symmetry, explicitly broken by quark masses; (ii) spon- 
taneous breakdown of chiral symmetry; (hi) analyticity, unitarity and crossing symmetry On 
the other hand, such a parametrization (effective Lagrangian) should be sufficiently general, 
and should not introduce any additional dynamical assumptions beyond those listed above, 
that could be hard to identify as emerging from QCD. The specificity of QCD then resides 
in numerical values of low-energy constants which characterize the above parametrization. 
The theoretical challenge is to calculate these low-energy parameters from the fundamental 
QCD Lagrangian. While such a calculation is awaited, these parameters can be subjected 
to experimental investigation. Chiral symmetry guarantees that the same parameters that 
are introduced through the low-energy expansion of QCD correlation functions also define 
the low-energy expansion of hadronic observables — pseudoscalar meson masses, transition 
and scattering amplitudes. 

In this paper, a new method will be elaborated that allows a detailed measurement 
of certain low-energy parameters, using the n — n elastic scattering data HH]. Instead 
of concentrating on a particular set of scattering lengths and effective ranges || whose 
extraction from experimental data is neither easy nor accurate, emphasis will be put on a 
detailed fit of the scattering amplitude in a whole low-energy domain of the Mandelstam 
plane, including the unphysical region. In this way it is possible to obtain some experimental 
insight on the low-energy parameter 2rhBo , where rh is the average of the up and down quark 
masses, B is the condensate 
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and |0) and F$ stand for the ground state and pion decay constant respectively at m u = 
ma = m s = 0. It is usually assumed that the parameter 2rhBo differs from the pion mass 
squared by not more than (1 — 2)% ||, and the standard chiral perturbation theory could 
hardly tolerate an important violation of this assumption On the other hand, this 
assumption has never been confronted with experiment otherwise than indirectly - through 
the Gell-Mann Okubo formula for pseudoscalar meson masses [[| . However, even the latter 
represents at best a consistency argument rather than a proof: the Gell-Mann Okubo formula 
can hold quite independently of the relation between 2mB and M% 0. An independent 
measurement of 2mBo is not only possible (as shown in the present work) but, for several 
reasons, it appears to be desirable: 

(i) The effective Lagrangian £ e ff contains, in principle, an infinite number of low-energy 
constants, which are all related to (gauge invariant) correlation functions of massless QCD. 
Among them, B plays a favored role: The order of magnitude of all low-energy constants 
other than Bq can be estimated using sum-rule techniques 0, which naturally bring in the 
scale A ~ 1 GeV characteristic of massive bound states. The expected order of magnitude 
of a low-energy constant related to a connected iV-point (N > 1) function of quark bilinears 
qYq, that is not suppressed by the Zweig rule or by a symmetry, is F$A 2 ~ N multiplied by 
a dimensionless constant of order 1. If quarks were not confined [Q, a similar estimate 
would relate B Q and the mass of asymptotic fermion states with quark quantum numbers. 
However, in a confining theory, no similar relation between Bq and the spectrum of massive 
bound states can be derived: qq is an irreducible color singlet and there is no complete set 
of intermediate states which could be inserted into the matrix element (0|gg|0). Bq could be 
as large as A ~ 1 GeV or as small as the fundamental order parameter of chiral symmetry 
breaking, F ~ 90 MeV. A priori, there is no way to decide in favor of one of these scales, 
at least before the non-perturbative sector of QCD is controlled analytically or by reliable 
numerical methods, using, for instance, sufficiently large lattices. In this paper, we suggest 
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how the question of the scale of B can be addressed experimentally. 

(ii) To the extent that 2rhB and (m s + m)B are close to M 2 and M^, respectively, the 
ratio of quark masses 

r = m s /m (1.2) 

M 2 

must approach — 1 = 25.9 There exists an independent measurement of the ratio 



r in terms of observed deviations from the Goldberger-Treiman relation (JTOj] in non-strange 
and strange baryon channels. This model-independent measurement indicates a considerably 
lower value for r than 25.9, unless the pion-nucleon coupling constant turns out to be below 



the value given by Koch and Pietarinen [TT[] by at least 4-5 standard deviations [[[2 

(hi) A reformulation of xPT which allows 2mS to be considerably lower than M 2 has 
been given in Ref. 0. It is as systematic and unambiguous as the standard xPT itself, and 
is particularly suitable in the case where Bq is as small as F n . It is based on a different 
expansion of the same effective Lagrangian, with the same infinity of independent terms. 
To all orders, the two perturbative schemes are identical but, in each finite order, they can 
(but need not) substantially differ. For each given order, the new scheme contains more 
parameters than the standard xPT, the latter being reproduced for special values of these 
additional parameters. Already at the leading order 0(p 2 ), the new scheme contains one 
additional free parameter 

2mB 

V = ^f (L3) 
If r] is set equal to 1, one recovers the leading 0(p 2 ) order of the standard %PT. The 

new expansion can therefore be formally viewed as a generalization of the standard scheme 

and - in this sense - it will be referred to as improved xPT, since it aims to improve the 

convergence of the standard perturbation theory. Demonstrating that such an improvement 

is irrelevant, by measuring, for instance, the ratio (1.3) and finding it close to unity, would 

be an important experimental argument in favor of the standard xPT. 

(iv) In some cases, the convergence of standard %PT actually appears to be rather slow. 

Most of the indications in this direction can be traced back to the fact that the leading 



0(p 2 ) order of the standard xPT underestimates the Goldstone boson interaction and, in 
particular, the n — n scattering amplitude. This manifests itself through virtual processes 
and/or final state interactions, as in 77 — ► 7r°7r° |I3"|| , 77 — 37r [[TJj], etc. It might even be 
that although the next order 0(p 4 ) improves the situation, it fails to reach the precision 
we may rightly expect from it. For example, the I — s wave n — it scattering length, 
which is <2q = 0.16 in leading order |I5| , gets shifted to a° = 0.20 by 0(p 4 ) corrections [0, 
while the "experimental value" |16HT^1 is a° = 0.26 ± 0.05. (In this paper it will be argued 



that scattering lengths are not the best quantities to look at. A more detailed amplitude 
analysis will reveal a possible amplification of the discrepancy, which exceeds one standard 
deviation.) 

(v) xPT should be merely viewed as a theoretical framework for a precise measure- 
ment of low-energy QCD correlation functions. Its predictive power rapidly decreases with 
increasing order in the chiral expansion : More new parameters enter at each order and 
more experimental data have to be included to pin them down. For this reason, a slow 
convergence rate might sometimes lead to a qualitatively wrong conclusion with respect to 
a measurement based only on the first few orders. This might concern, in particular, the 
measurement of the ratio 77 (1.3) within the standard xPT. In the corresponding leading 
order, 77 is fixed to be 1, independently of any experimental data. This property of standard 
%PT could bias the measurement of 77 if 77 turned out to be considerably different from 1: 
one would presumably have to go to a rather high order and include a large set of data to 
discover the truth. In this case, the improved xPT would be a more suitable framework to 
measure 77 faithfully. The reason is that in the improved xPT, 77 is a free parameter from 
the start: It defines the leading order tc — tt amplitude. Neglecting, for simplicity, Zweig-rule 
violation (cf. Ref. and Sec. IV A), the latter reads 

A(s\tu) = -^(s- V M 2 ). (1.4) 

Using in this formula the value of a° = 0.26 ± 0.05, one concludes that 77 = 0.4 ± 0.4 already 
at the leading order. The measurement then has more chances to saturate rapidly - say, 
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at the one-loop level - provided rj is much closer to 0.4 ± 0.4 than to 1. The same remark 
applies to measurements of the quark mass ratio r = m s /m, which, incidentally, is closely 
related to rj |K|: A slow convergence of the standard xPT could lower the leading-order 
result r = 25.9 by considerably more than the usually quoted (10 — 20)% P,|19 . 



(vi) The question of the actual value of r] and/or of r = m s /rh has to be settled ex- 
perimentally. None of the known properties of QCD, nor the fact that light quark masses 
are tiny compared with the hadronic scale A ~ 1 GeV, imply that rj should be close to 1 
and that r should be close to 25.9. The proof of this negative statement is provided by the 
existence of a mathematically consistent generalization of the standard xPT that does not 
contradict any known fundamental property of QCD and allows for any value of rj between 
and 1 (and for any value of r between 6.3 and 25.9) ||. Only in the special case of r\ and r 
close to 1 and 25.9, respectively, can the standard %PT claim a decent rate of convergence. 

In Sec. II, the precise mathematical definition of the improved xPT, in terms of the 
effective Lagrangian, is briefly summarized. It is not a purpose of this paper to present a 
full formal development of this theory; incidentally, most of it can be read off from existing 
calculations || after rather minor extensions (which will be presented elsewhere). Here, 
we will mainly concentrate on phenomenological aspects of the problem in connection with 
low-energy tt — tt scattering. The content of Sec. Ill is independent of any particular 
xPT scheme. In that section, a new low-energy representation of the tt — tt scattering 
amplitude is given that provides the most general solution of analyticity, crossing symmetry 
and unitarity up to and including the chiral order 0(p 6 ). (Partial wave projections of this 
representation coincide with a particular truncation of the well-known Roy equations @.) 
Subsequently, this representation is used both to constrain the experimental data and to 
perform a comparison with theoretical amplitudes as predicted by the two versions of x?T. 
For the case of the improved %PT, the one-loop amplitude is worked out in Sec. IV. Finally, 
a method permitting a detailed fit of the experimental amplitude in a whole low-energy 
domain of the Mandelstam plane is developed in Sec. V. This method is then applied to 
various sets of existing data. 



II. FORMULATION OF IMPROVED xPT 



Following rather closely the off-shell formalism which was elaborated some time ago by 
Gasser and Leutwyler |3|], we consider the generating functional Z(v^,a^,x) °f connected 
Green functions made up from SU(3) x SU(3) vector and axial currents as well as from 
scalar and pseudoscalar quark densities, as defined in QCD with three massless flavors. The 
sources v M , a M and x are specified through the Lagrangian 

£ = Cqcd + ^7s)g - qnXQL - QlX^QRi (2.1) 

which defines the vacuum-to- vacuum amplitude expiZ. Here, q^R = |(lT75)s stand for the 
light quark fields u, d, s and Cqcd is invariant under global SU(3) x SU(3) transformations 
of qi and q^. and are traceless and hermitean, whereas 

X = s + ip (2.2) 

is a general 3x3 complex matrix (s and p are hermitean). Explicit chiral symmetry breaking 
by quark masses is accounted for by expanding Z around the point 

/ \ 



v » = a** = 0, X = M q = 



(2-3) 



y m 8 J 

The scalar-pseudoscalar source x an d the quark mass matrix M. q are closely tied together 
by chiral symmetry. (Notice that our source x differs from the x defined in Ref. [[J by a 
factor of 2B .) 

Instead of calculating Z, the effective theory parametrizes it by means of an effective 
Lagrangian which depends on the sources and on eight Goldstone boson fields 

(/(x)=exp^%(4 (2.4) 
0= i 

Leaving aside anomaly contributions described by the Wess-Zumino action, the effective 
Lagrangian C e ff(U, v M , a M , x) is merely restricted by the usual space-time symmetries and 
by the requirement of invariance under local chiral transformations [£Il,r G SU (3)] 
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U(x) - Q R (x)U(x)Ql(x), x (x) - n R (x) X (x)Ql(x) (2.5) 
compensated by the inhomogeneous transformation of the sources and a M : 

u" - a' 1 -> - + id")^. (2.6) 

(This gauge invariance of the nonanomalous part of Z is necessary and sufficient to re- 
produce all SU(3) x SU(3) Ward identities.) Otherwise, the effective Lagrangian remains 
unrestricted. 

C e ff can be written as an infinite series of local terms, 

£eff = ^ ^ nm £nm, (2.7) 
n,m 

where C nm denotes an invariant under the transformations (2.5) and (2.6) that contains the 
n-th power of the covariant derivatives and the m-th power of the scalar-pseudoscalar 
source x- The sum over independent invariants that belong to the same pair of indices (n, m) 
is understood. The covariant derivatives are defined as 

D„U = d^U - i(vp + a^U + iU( Vlx - O, (2.8) 

and likewise for D^x- The expansion coefficients £ nm represent properly subtracted linear 
combinations of massless QCD correlation functions that involve n vector and/or axial cur- 
rents and m scalar and/or pseudoscalar densities, all taken at vanishing external momenta. 
The first two terms in the sum (2.7), for instance, read (n is even) 

f°C 20 = - A F%{D^D,U). (2.9) 

Everything said so far is rather general and independent of any particular perturbative 
scheme. 

Chiral perturbation theory is an attempt to reorder the infinite sum (2.7) as 
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£e// = E£ W > (2-10) 

d 

where collects all terms that in the limit 

p^O, Af ff -»•(), p 2 /Ml fixed (2.11) 

behave as 0{p d ) (p stands for external momenta). In order to relate the expansions (2.10) 
and (2.7), one needs to know the effective infrared dimension d{m q ) of the quark mass. The 
invariant C nm then contributes as 0(p dnm ), where 

d nm = n + m d(m q ) . (2.12) 

For infinitesimally small quark masses, one should have 

d(m q ) = 2, m q -> 0. (2.13) 
This follows from the mathematical fact that in QCD 

lim K \^ )B ° = 1. (2.14) 



(Here, i,j — u,d,s and M P is the mass of the pseudoscalar meson Tj, i ^ j.) The assumption 
that in the real world, i.e., for physical values of quark masses, the effective dimension of 
the quark mass is 2, underlies the standard xPT. It amounts to the well-known rule which 
asserts that each insertion of the quark mass matrix and/or of the scalar-pseudoscalar source 
X counts as two powers of external momenta. Equivalently, the standard can be viewed 
as an expansion around the limit 

(p,m ff )->0, p 2 /M 2 P fixed. (2.15) 

Since, by definition, the low-energy constants £ nm are independent of quark masses, they are 
0(1) in the limit (2.15). 

It is easy to see that the convergence of the standard xPT could be seriously disturbed if 
Bo€A~l GeV, say B ~ 100 MeV [p|,PH|. The expansion of Mp reads (i,j — u,d,s; i ^ 

j) 
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Mp = (m, + mj)B + (m* + m/Ao + . . . , (2.16) 

where the dots stand for non-analytic terms and for higher order terms. Aq can be expressed 
in terms of two-point functions of scalar and pseudoscalar quark densities divided by Fq 
TUfl . It satisfies a superconvergent dispersion relation, whose saturation leads to the order 



of magnitude estimate A ~ 1 — 5. For B as small as 100 MeV, the first and second order 
terms in Eq. (2.16) then become comparable for quark masses as small as (10 - 50) MeV. 
In order to accommodate this possibility, the improved x_PT attributes to the quark mass 
and to the vacuum condensate parameter B the effective dimension 1, 

d(m q ) = d(B ) = 1, (2.17) 

reflecting their smallness compared to the scale A. This does not contradict mathematical 
statements such as (2.14). It only means that, although for physical values of quark masses 
the ratio in Eq. (2.14) remains on the order of 1, it is allowed to differ from 1 considerably. 

To summarize, in the improved %PT each insertion of the quark mass-matrix Ai q and/or 
of the scalar-pseudoscalar source \ counts as a single power of external momentum (pion 
mass) and so does the parameter B . This leads to a new expansion of the effective La- 
grangian 

C eff = J2C d , (2.18) 

d 

where each C d contains more terms C nm than does the corresponding term C d in the case of 
the standard counting. The improved xPT is a simultaneous expansion in p/A,m q /A and 
Bq/A around the limit 

{p,m q ,B )^0, p 2 /M 2 P and m q B /M P fixed. (2.19) 

This is just another way to realize the chiral limit (2.11). The fact that - in the effective 
theory - we treat Bq as an arbitrary expansion parameter does not contradict the general 
belief that, within QCD, this parameter is fixed and - hopefully - calculable. After all, 
quantum electrodynamics is also based on an expansion in a, in spite of the general belief 
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that there might exist a more fundamental theory in which the value of a is fixed and 
calculable pi ]. 



III. RECONSTRUCTION OF THE LOW-ENERGY 
it -ir SCATTERING AMPLITUDE 
NEGLECTING 0(p 8 ) EFFECTS. 

The analysis of low-energy tc — tc scattering, traditionally based on analyticity, crossing 



symmetry and unitarity |22|,|T6|,|20| , considerably simplifies if, in addition, one takes into 
account the Goldstone character of the pion. First, in the chiral limit (2.11), higher (£ > 
2) partial waves are suppressed. The reason stems from the fact that in the limit (2.11) 
the whole amplitude behaves as 0(p 2 ) and furthermore, it does not contain light dipion 
bound state poles. Unitarity then implies that the scattering amplitude is dominantly real, 
since its imaginary part behaves as 0(p 4 ). Analyticity then forces the leading 0{p 2 ) part 
of the amplitude A(s\tu) to be a polynomial in the Mandelstam variables. Furthermore, 
higher than first order polynomials are excluded: They would be 0(p 2 ) only provided their 
coefficients blew up as M 2 — > 0, which would contradict the finiteness of the S-matrix in 
the limit m q — > with the external momenta kept fixed at a non-exceptional value. Finally, 
crossing symmetry allows one to express the 0{p 2 ) part of the scattering amplitude A(s\tu) 
as 

A Lead (s\tu) = ^- 2 M 2 + ^- 2 (3s - AM 2 ), (3.1) 

where a, f3 are two dimensionless constants which are 0(1) in the chiral limit. The linear 
amplitude (3.1) does not contribute to i > 2 partial waves. Consequently, the latter behave 
in the chiral limit as 0(p 4 ) and, owing to unitarity, the absorptive parts of £ > 2 waves 
are suppressed at least to 0(p 8 ). This conclusion holds independently of more quantitative 
predictions of x?T, which in the actual case merely concern the values of the two parameters 
a and (3 in Eq. (3.1). 
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The second simplification resides in the suppression of inelasticities arising from inter- 
mediate states that consist of more than two Goldstone bosons. The behavior of ra-pion 
invariant phase space in the chiral limit (2.11) is given by its dimension: It scales like p 2n ~ 4 . 
Amplitudes with an arbitrary number of external pion legs are dominantly 0(p 2 ). Conse- 
quently, the contribution of multi-pion (n > 2) intermediate states to the absorptive part of 
the elastic tt — n amplitude is suppressed in the chiral limit at least to 0(p s ). 

The smallness of higher partial waves and of inelasticities are of course well-known phe- 



nomenological facts [16]. It is important that these "remarkable accidents" (see page 53 of 
Fj| ) can be put under the rigorous control of chiral power counting: The previous discussion 
suggests that a rather simple amplitude analysis of low-energy it — it scattering can be per- 
formed up to and including 0(p 6 ) contributions. In the following we confirm and elaborate 
this expectation in detail. It will be shown in particular that, neglecting 0(p s ) contributions, 
the whole scattering amplitude can be expressed in terms of low-energy s and p wave phase 
shifts and six (subtraction) constants. (The latter are related to the experimental phase 
shifts via unitarity.) The resulting expression (3.2) will prove particularly useful both for 
constraining low-energy experimental data and for providing a basis for a confrontation of 
chiral perturbation theory up to two loops with experiment. 

A. Statement of the theorem 

Let A denote a scale (slightly) below the threshold for production of non-Goldstone 
particles. The 7r — ir amplitude can be written as 

3 -A(s\tu) = T(s) + T(t) + T(u) 



327T 

+ l[2U(s)-U(t)-U(u)} 

+ ±[( s -t)V{u) + {s-u)V{t)} 

+ R A (s\t,u). (3.2) 
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The remainder, R\, behaves in the chiral limit as 0(p 8 ) relative to the scale A: up to possible 
logarithmic terms, 

R A = 0([p/A} 8 ), (3.3) 

where p stands for external pion momenta. In practice, A < 1 GeV. The functions T, U 
and V are analytic for s < iM^., whereas for AM 2 < s < A 2 their discontinuities are given 
by the three lowest partial wave amplitudes 

ImT(s) = l -{Imf (s) + 2Imf (s)} 

ImU(s) = i{2/m/ °( S ) -5Jm/ 2 (5)} 
27 1 

ImV(s) = - —^ Imft(s). (3.4) 

The real parts of the functions T, U and V are defined only up to polynomials 

6T(s) = x(s - ^Ml) 

SU (s) =y + yis + y 2 s 2 + y 3 s 3 

5V(s) = -{y t + AMly 2 + IQM%) + (y 2 + UM*y 3 )s - 3y 3 s 2 , (3.5) 

where x and the y's are five arbitrary real constants: because of the relation s + t + u = 4M 2 , 
the two sets of amplitudes T, U, V and T + 5T, U + 8U, V + 5V lead to the same scattering 
amplitude A. (It is shown in Appendix B that Eqs. (3.5) actually represent the most general 
transformation of T, U, V leaving the scattering amplitude invariant.) After conveniently 
fixing the "gauge freedom" (3.5), the functions T, U and V can be written as 

T(s) =t + t 2 s 2 + t 3 s 3 + — [ A ^^—ImT(x) 

7T J4A/2 X 6 X — S 

U(s)= - ^—MU(x) 

71 J AMI X A X — S 

s 2 r^ 2 dx 1 

V(s) = v 1 + v 2 s + v 3 s 2 + — / — ImV(x), (3.6) 

7T HMl X z X — S 



Notation and normalization are reviewed in Appendix A. 
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where the imaginary parts are given by Eqs. (3.4) and the t's and i>'s are constants. It 
will be shown shortly that Eq. (3.2) is a rigorous consequence of analyticity and crossing 
symmetry and of the Goldstone nature of the pion. 



B. Unitarity 

The low-energy representation (3.2) of the scattering amplitude is exact up to an 0(p 8 ) 
remainder. In the whole interval 4M 2 < s < A 2 , unitarity can be imposed with the same 
accuracy in terms of partial waves //. As already pointed out, deviations from the unitarity 
condition 



above the inelastic threshold are of the order 0(p 8 ). The amplitude (3.2) contains all partial 
waves. For £ > 2, the partial waves are real. Nevertheless, unitarity automatically is 
satisfied for £ > 2 up to 0(p 8 ) terms, since higher partial waves anyway are 0(p 4 ) or smaller. 
Consequently, it is sufficient to impose unitarity for the three lowest waves /° , fl and /q 
(hereafter denoted as f a , a = 0, 1, 2 according to their isospin). Projections of Eq. (3.2) into 
the three lowest partial waves read 

7T J 4MJ X 6 X — S 
"A 2 dx 2 



1 r A dr 1 

+ - W ^ x)Im f b (x) + 0(p 8 ) (3.7a) 

7T UmI X b=Q 

for the two s waves (a = 0, 2), whereas the p wave projection is 

p f n dm I * 2 (*-4M?) f A2 dx Imh(x) 
Re}i(s) = Pi(s) H f — 

7T J iM%X 2 (x - AMI) X ~ S 

+ - / — E A(^) + 0(/) (3.7b) 

Here, -P a (s) are third order polynomials whose coefficients are defined in terms of the six 
constants to, t 2 , t 3 and vi, v 2 , v 3 which appear in Eqs. (3.6). These polynomials are tabulated 
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in Appendix C, together with the nine kernels W a b(s,x) which define the left-hand cut 
contributions to the partial waves. 

Eqs. (3.7a) and (3.7b) may be viewed as a particular truncation of the infinite system of 
Roy equations, which slightly differs from the form in which these equations have been used 
in the past [23|. Here, the truncation in angular momentum and energy is performed under 
the systematic control of chiral power counting. In particular, Eqs. (3.7a,b) do not require 
a mo del- dependent evaluation of "driving-terms" which in the standard treatment behave 
in the chiral limit as 0(p 4 ), owing to the use of twice-subtracted dispersion relations. The 
price to pay is the occurrence of six (a priori unknown) constants in the polynomials P a (s) 
instead of only two constants (usually, the two s wave scattering lengths) which characterize 



the inhomogeneous terms in standard Roy equations f2(j|23| . 

Eqs. (3.7a) and (3.7b) can be used to fully reconstruct from the data the whole am- 
plitude A(s\tu) up to and including accuracy 0(p 6 ) in the whole low-energy domain of the 
Mandelstam plane, including the unphysical region. For this purpose one has to know the 
absorptive parts of three lowest partial waves for 4M 2 < s < A 2 and the six constants 
t , t 2 , t 3 , V2, v%, v 3 . Suppose one knew Im f a (s) with associated error bars in the whole inter- 
val 4M 2 < s < A 2 < 1 GeV 2 . Then one could calculate the dispersion integrals on the right 
hand side of Eqs. (3.7a,b). One would then determine the constants t and v from the best fit 
to the values Re f a (s) determined from the input Im f a (s) via the unitarity condition. The 
X 2 of this fit may be considered as a measure of the internal consistency of the input data 
Imf a (s). In practice, experimental information on Im f a (s) is only available for s well above 



the threshold. In this case, a more sophisticated iteration procedure plfof Eqs. (3.7a,b) 
has to be used in order (i) to extrapolate the experimental data down to the threshold and, 
simultaneously, (ii) to determine the six constants t and v. In both cases, the resulting 
amplitude is given by the formula (3.2). 
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C. Proof of the reconstruction theorem 



Formulae (3.2) and (3.6) can be proven following the original derivation of the Roy 
equations p0| . The proof is based on fixed t dispersion relations for the three s-channel 
isospin amplitudes 



F(s,t,u) 



(s,t,u), (3.8) 



combined with the crossing symmetry relations 

F(s, t, u) = C su F{u, t, s) = C st F(t, s, u) = C ut F(s, u, t). (3.9) 

(Properties of the crossing matrices C su , C s t and C u t are reviewed in Appendix A.) The 
standard Roy equations are derived from twice-subtracted dispersion relations - cf. the 
minimal number of subtractions required by the Froissart bound. In this case, however, the 
high-energy tail of the dispersion integral, which is hard to control in a model independent 
way, contributes to the 0(p 4 ) part of the amplitude. (In standard Roy equations, this 
contribution is contained in the so-called driving terms ||23|| .) If, on the other hand, one 
requires at low energy the precision 0(p 4 ) or higher, then it is more appropriate to stick to 
less predictive triply- subtracted fixed-t dispersion relations: 

F(s, t) = C st {8i + (t) + (s- it)b_(t) + (s- u) 2 c + (t)} 

+ - r + — C s i/mF(M). (3.10) 



7T J4A/| X 3 [X — S X — U 

Here the subscript ± refers to the eigenvalues ±1 of the crossing matrix C tu - [Notice that in 
the s-channel isospin basis (3.8), Ct u = diag(+l, —1, +1).] The subtraction term then rep- 
resents the most general quadratic function in s (for fixed t) symmetric under s — u crossing. 
By construction, the dispersion integral in Eq. (3.10) exhibits s — u crossing symmetry too. 
The task is now to impose the remaining two crossing relations and to determine the sub- 
traction functions a,b, and c. This can be achieved, neglecting in Eq. (3.10) contributions 
of chiral order 0(p 8 ) and higher. 

17 



Let A be a scale set by the threshold of production of non-Goldstone particles. Let us 
split the dispersion integral in Eq. (3.10) into low energy (x < A 2 ) and high energy (x > A 2 ) 
parts. For 4M 2 < s < A 2 , the imaginary part can be written as 



ImF = Im& + (s) + 1 + 



2t 



7m$_(s) + A e > 2 (s,t), 



s - 4M2 J 

where the first two terms stand for the contributions of s and p waves: 

/ „ \ 



(3.11) 



Im& + (s) 



V 



/ 







3Imfl(s) 




(3.12) 





Imflis) 

A^> 2 then collects the absorptive parts of all higher partial waves. The reason for this 
particular splitting resides in the chiral counting mentioned at the beginning of this section: 
The first two terms in Eq. (3.11) dominantly behave as 0(p 4 ), whereas A^> 2 is suppressed 
to 0(p 8 ). The dispersion integral I(s, t) in Eq. (3.10) then splits into three parts, 



I(s, t) = I e<2 (s, t) + Lj> 2 (s, t) + I H (s, t). 



(3.13) 



I^<2 (\t>2) is the contribution of low-energy t < 2 (£ > 2) partial waves, and I# represents 
the high frequency part in which no partial wave decomposition is performed. Extracting 
from 1h its leading low energy behavior, one can write 



I H (s, t) = (s 3 + u 3 C su )B. A + R H , 



(3.14) 



where Ha are constants which can be expressed as integrals over high-energy n — n total 
cross sections, and the remainder behaves at low energies as 



R„ = 0([p/A] 8 ). 



(3.15) 



The low-energy high angular momentum part Ig>2 is also suppressed to 0(p 8 ), reflecting the 
leading behavior of the absorptive part A^> 2 in the chiral limit and the fact that the cor- 
responding dispersion integral (3.10) extends over a finite interval x G [4M 2 , A 2 ]. Hence, it 
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remains to concentrate on the low-energy low angular momentum part le<2- Using Eq. (3.11), 
one easily checks the identity 



and <l>± denote the following dispersion integrals over the imaginary parts of low-energy s 
and p waves [cf. Eq. (3.12)]: 



One observes from Eq. (3.17) that the function <f>(s,t,u) exhibits the full three-channel 
crossing symmetry. Furthermore, the second and third terms in Eq. (3.16) represent a 
function that is quadratic in s (at fixed t) and symmetric under s — u crossing. These terms 
can therefore be absorbed into the subtraction polynomial in the dispersion relations (3.10) 
by a suitable redefinition of (yet unknown) subtraction functions a + , b , c+. Consequently, 
the whole amplitude F can be rewritten as 




(3.16) 



where 




(3.17) 




(3.18) 



F(s, f) = *(s, t, u) + P(s, t, u) + 0([p/A] 8 ), 



(3.19) 



where P is of the form 



P = C st {a+(t) + (s- u)f3_(t) + (s- uf 1+ {t)} 
+ (s 3 + m 3 C su )H a . 



(3.20) 
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Notice that the unspecified 0(p s ) contributions in Eq. (3.19) originate both from the high- 
energy remainder R# (3.14) and from the low-energy higher angular momentum part Lj>2- 
Crossing symmetry of the scattering amplitude F should hold order by order in the chiral 
expansion. Since the function $ (3.17) exhibits full crossing symmetry, it remains to impose 
the latter for the function P (3.20). Because of the manifest s — u symmetry, it is enough 
to require 

P(s,t,u) = C st P(t,s,u). (3.21) 

Neglecting 0(p 8 ) contributions, this equation represents the necessary and sufficient condi- 
tion for the complete crossing symmetry of the amplitude F. 

Eq. (3.21) can be easily solved. Considering s and t as independent variables, one easily 
finds that a+(t), (3-(t), and j+(t) should be cubic, quadratic and linear functions of t 
respectively. Hence, P(s, t, u) is a general crossing symmetric polynomial in the Mandelstam 
variables of (at most) third order. Such a polynomial contains six independent parameters 
(see Appendix A). Indeed, after some simple but lengthy algebra, one verifies that Eq. (3.21) 
leaves a six parameter freedom in the original expression (3.20) for P. 

It remains to rewrite the result (3.19) in terms of the single amplitude 

A(s\tu) = A(s\ut) = ^ {f (0 \s, t, u) - F {2 \s, *,«)}. (3.22) 

The function <fr gives rise to a contribution of the form (3.2) in which only the disper- 
sion integrals of Eq. (3.6) occur. (One easily checks that ImT, ImU and ImV are given 
by Eqs. (3.4).) Furthermore, taking into account the ambiguity (3.5) in the definition of 
T, U and V, it is clear that a general crossing symmetric polynomial may be conveniently 
parametrized by the six independent parameters t ,t 2 , t 3 , i>i, v 2 , v 3 as in Eqs. (3.6). 

IV. PERTURB ATIVE tt - tt AMPLITUDE AND THE 
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EFFECTIVE INFRARED DIMENSION OF THE QUARK MASS 



We are now in a position to compare the two alternative low-energy expansions of the 
amplitude A(s\tu) generated by chiral perturbation theory according to the two possible 
values of the effective dimension of the quark mass: 2, in the case of the standard x?T, 
and 1 in the case which was defined in Sec. II as improved x?T. Up to and including two 
loops, the amplitude A should be of the general form (3.2). Consequently, neglecting 0(p s ) 
contributions, one can work with the three functions T, U and V of a single variable and 
decompose them as 

T(s) = £ TW( a ), U(s) = £ <»)(*), V(s) = £ V^(s), (4.1) 

n=0 n=0 n=0 

where n refers to the number of loops (including tree contributions of the corresponding 
order). It will be shown that the amplitudes T, U and V start to be sensitive to the effective 
dimension of the quark mass at leading (n = 0), one-loop (n = 1) and two-loop levels 
respectively. 



A. Leading 0(p 2 ) order 

If the dimension of the quark mass is 2, i.e., if each power of the scalar pseudoscalar source 
X in £ e ff counts for two powers of pion momentum (mass), then the effective Lagrangian is 
dominated by the well-known expression 

= ±F*{{(D«U) + (D„U)) + 2B (x + U + U + X )}- (4.2) 

This formula collects all possible invariants of dimension 2. To leading order, the pion and 
the kaon masses read 

Ml = 2mB 

M\ = (m s + m)B (4.3) 



and the ir — n amplitude takes the well-known form, first given by Weinberg |T5| 
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A lead (s\tu) = -2(3 - 2mB ) = -^(s - M 2 n ) (4.4) 

This represents the standard scenario of chiral perturbation theory. It can hardly be circum- 
vented provided the scale of B is large compared to the pion mass, typically, B > 1 GeV. 

On the other hand, if B turned out to be much smaller than the GeV-scale, e.g., 
comparable to the fundamental order parameter F (~ 93 MeV), then the above way 
of counting effective infrared dimensions would be modified. Both the quark mass and 
the condensate B should then be considered as quantities comparable to the pion mass. 
They should both be attributed effective infrared dimension 1 and they should both be 
viewed as expansion parameters. In this case, every insertion of the source x( x ) counts as a 
single power of pion momentum and the formula (4.2) no longer represents the most general 
expression of dimension 2. Instead, the complete collection of invariants of dimension 2 now 
reads 

= - A F*{{D»UD^U + ) + 2B ( X + U + X U + ) 

+ A (x + U X + U + x u + xu + ) + z s Q { x + u + u + x ) 2 + 

+ Z£(x + U - X U + ) 2 + 2H (x + x)}- (4.5) 

where the tilde over the symbol C here (and below) indicates the use of the modified chiral 
power counting. The terms containing two powers of x are usually included into the next-to- 
the-leading part of the effective Lagrangian. Here, they appear of the same dimension 
and they are expected to be of a comparable size as the standard expression (4.2). The 
low-energy constants A , Zq and Zq represent appropriately subtracted zero-momentum 
transfer two-point functions of scalar and pseudoscalar quark densities, divided by Fq. These 
two-point functions are order parameters of spontaneous chiral symmetry breaking and, 
consequently, they satisfy super convergent dispersion relations. A simple saturation of the 
latter with a few of the lowest massive hadronic states suggests that the dimensionless 
constants Aq and Zq are of the order 1, say, Aq ~ 1 — 5. On the other hand, Z violates 
the Zweig rule in the ++ channel and consequently it is expected to be suppressed. The 
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parameters Zq , Zq and Aq are related to the low-energy constants Lq, L-j and Lg of the 
standard d = 4 Lagrangian ||. Expanding the latter constants in powers of B , one 
getsQ 

L«=^£){Z% + { Bl)l 




Ls=(-§^) 2 {Ao + 0(Bt)}, (4.6) 

where O(Bq) terms represent divergent contributions to the two-point functions defining the 
divergent parts of the bare constants L 6 , L 8 . (The constants A , Zq and Zq do not undergo 
any infinite renormalization.) 

The leading order pion and kaon masses (denoted by a tilde) now read 

Ml = 2m{B + 4mZ 5 ) + 4m 2 A 

M\ = (m s + m)(B + AitiZq) + (m s + m) 2 A . (4.7) 

Here B stands for the dominant 0(p) contribution to the SU(2) x SU(2) quark- ant iquark 
condensate (divided by Fq) taken at m u = rrid = 0: 

(uu) mu=md=0 = (dd) mu=md=0 = -FqB + 0(m 2 s ). (4.8) 

Within the modified chiral power counting, B consists of two terms 

B = B + 2m s Z^ (4.9) 

which are both of the order 0(p). In principle, they could be of comparable size, if Zq were 
not suppressed by the Zweig rule. 

2 The order of magnitude estimate Aq ~ 1 — 5 is compatible with the standard xPT estimates. 
Taking Aq ~ 5, and using the standard value Bq ~ 1.2 GeV, the ^-contribution to Lg in Eq. (4.6) 
becomes 1.6 x 1CP 3 , which is consistent with the standard xPT measurement of Lg p|. 
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The leading contribution to the n — it scattering amplitude calculated from the improved 
0(p 2 ) Lagrangian (4.5) turns out to be independent of low-energy parameters A and Z$ , 
and it can be expressed in term of the quark- ant iquark condensate B, 

A lead {s\tu) = i_( s -2m£), (4.10) 

in complete analogy with the standard result (4.4). Although Eq. (4.10) and Weinberg's 
formula (4.4) formally coincide if one neglects Zweig-rule violation, their numerical content 
is rather different, because of different scales of quark- ant iquark condensation in each %PT 
alternative. In Eq. (4.4), 2mBo is the leading approximation to M%, whereas in the improved 
%PT, the relation between the quark-antiquark condensate and the pion mass is more subtle: 
Indeed, using first Eq. (4.7), formula (4.10) can be rewritten as 

1 _ m 2 

A lead (s\tu) = -5 (s - Ml) + -f 6 (1 + 2C), (4.11) 

where 

4rh 2 A Zl 



C=-f. (4.12) 



Ml ' A 

Whereas in the standard x?T e would be a small quantity of the order 0(p 2 ), in the improved 
%PT, e is 0(1) and there is no reason for it to be particularly small; hence, the second term 
in Eq. (4.11) represents a leading order modification of the Weinberg's formula (4.4). (£ 
measures the Zweig rule violation in the ++ channel and can be expected rather small.) 
Using Eqs. (4.7) one may easily check that e can indeed be of order 1 for natural values of 
Aq (cf. footnote |2]) and for reasonably small values of quark masses. Setting — for the sake 
of illustration — Bq = 150 MeV and rh = 25 MeV, and neglecting Zweig rule violation, 
one obtains e = 0.62, A = 4.8 and m s ~ 195 MeV . 

The leading order mass formula (4.7) implies a relation between the parameter e and the 
quark-mass ratio r = m s /m: 



r 2 — r M j 



r 2 = 2^- - 1 ~ 25.9 (4.13) 



r z - 1 M 2 
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If r decreases from its canonical leading order value r = r2, then e increases and reaches 1 
for r = r 1; 

n = 2^fc - 1 ~ 6.33. (4.14) 

M n 

Similarly, the order parameter B can be expressed as 

.a = 1 _ [1 + ( r + 2)C]e. (4.15) 

7T 

This ratio decreases from its canonical value 1 down to zero, as r decreases from r = r 2 to 
f — r crit{() ~ ri, for which i?o vanishes. Notice that stability of the massless QCD vacuum 
under perturbation by small quark masses implies B > 0. 

B. Next to the leading 0(p 3 ) contribution 

In the improved chiral perturbation theory, the leading order Lagrangian is followed 
by a dimension 3 term C^ 3 \ which contributes at the tree level before one-loop contributions 
of dimension 4 start to appear. £ {3) reads 

^ (3) = \f^{d^u + d^ x + d,x + d^u) 

+ pi{{x + uf + ( x u + ) 3 ) + P 2(x + x(x + u + u + x )) 

+ P3((X + U) 2 - (xU + ) 2 )(x + U - X U + ) + ...}. (4.16) 

The dots stand for terms that violate the Zweig rule in a nonanomalous channel. Notice that 
(4.16) differs in its first term from the expression given for in Ref. |J. The two forms 
of are equivalent: they are related by a simple redefinition of the Goldstone boson field 
U. The low energy constants £ and pi are finite — there are no divergences of dimension 3. 
induces a shift in the pion mass, 

5Ml = eM^(9A! + A 2 ), (4.17) 

where 
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* = (4,8, 

are dimensionless parameters of order 0(M n ). Similarly, the leading tt — tc amplitude receives 
a constant d = 3 contribution 

5A(s\tu) = (81 Ax + A 2 ). (4.19) 

Finally, the first term in is responsible for splitting of the decay constants F n , Fx, F^. 
Eliminating the low-energy parameter £, one obtains, to that order 

^ = 1 + -!)■ ( 4 -20) 

It is convenient to collect all d = 2 and d = 3 contributions, and to express the resulting 
tree amplitude in the form (3.1): 

1 M 2 
A tree (s\tu) = -^[ocMl + f3(3s - AMI)} + ^f 2 8a, (4.21) 

where M n and F w denote the experimental (charged) pion mass and decay constant The 
parameters a and (3 read 

^ = 1 + 36(1 + 20, 0=S (4.22) 
P -^o 

whereas 5a = £0:3 + 5«4 describes small 0(p 3 ) and 0(p 4 ) corrections. 5a arises from the 

genuine 0(p 3 ) contribution (4.19) of C^ 3 ' to the it — it amplitude and from the introduction 

of the physical mass into the formula (4.21). Using Eqs. (4.17) and (4.19), the 0{p 3 ) 

constant M% 5a 3 can be expressed in terms of the parameters Ai and A2 of C^ 3 ': 

Ml 6a 3 = e(3Ml [72Ai - {27\ x + 3A 2 ) e (1 + 2C)]. (4.23) 

The remaining term M 2 5a^ accounts for the 0(p 4 ) and higher contributions to M 2 . One 
has 



3 In practice, = 139.6 MeV and F w = 93.1 MeV will be identified with the corresponding 
theoretical expressions up to and including the highest order of xPT considered. 
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M^a 4 = -/3AA^[l + 3e(l + 2C)], (4.24) 

where 

AM^ = Ml-Ml- 5Ml (4.25) 

represents the 0(p A ) difference between the physical value and the tree approximation of 
the pion mass squared. 

The results of standard are reproduced by setting e = ( = in the previous 

, o 

equations; i.e., r = r 2 ~ 25.9. In this case, M% reduces to M\ [Eq. (4.3)], and a — (3 ~ 1. 
The improved %PT still requires f3 ~ 1, but a is now allowed and expected to be considerably 
larger, since e is now an 0(1) quantity. In fact, the vacuum stability conditions mentioned 
above imply that for a given quark mass ratio r [lying between r\ and r 2 - cf. Eqs. (4.13) 
and (4.14)], the Zweig rule violating parameter ( = Zq /A should satisfy 

Using these bounds in Eq. (4.22), one obtains a rather narrow band of allowed values in the 
plane defined by the ratio a/ (3 and r. This band is shown in Fig. I. 

It is straightforward to rewrite the above result in terms of the amplitudes T, U and V. 
The tree contribution to these amplitudes simply reads 

r(°'(s) = (& + 5a)Ml, U i0 \s) = 0, \/ {0) (s) = 9/3, (4.27) 

where 

and likewise for 5a. Our main task is to use all available experimental information to 
measure a, f3 and, indirectly, the quark mass ratio r. 

C. One loop 0(p A ) order 



Let C nm denote an invariant entering the effective Lagrangian, that contains n powers of 
covariant derivatives D and m insertions of the scalar-pseudoscalar source x- (F° r simplicity, 
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the expansion coefficients £ nm of Eq. (2.7) are included in C nm .) In the standard chiral 
perturbation theory the d — 4 part of the effective Lagrangian can be written as 

£ (4) = £ £nm- (4.29) 

n+2m=4 

It contains all counterterms which are needed to renormalize one-loop contributions gener- 
ated by C^ 2 '. If the dimension of the quark mass is 1, one-loop renormalization gets modified 
in two respects: (i) The effective dimension of a term L nm is d = n+m instead of d = n+2m, 
and (ii) Bq is now a (small) expansion parameter of dimension 1. It follows, in particular, 
that renormalization has to be performed order by order in B . The modified d = 4 part of 
C e ff then reads 

£ (4) = J2 £n m + B (£ 21 + £ 03 )+B 2 £ 02 . (4.30) 

n+m=4 

The last two counterterms are needed to renormalize the -Bo-dependent part of one-loop 
divergences generated by C^ 2 \ Terms which are contained both in and in are 
merely made with four derivatives 

£ 40 = L^D^U+D^U) 2 + L 2 (D^U + D U U){D»U + D U U) 
+ L Z {D^U + D^UD U U + D V U) 
- iL 9 (F*D f *UD u U + + F^D»U + D U U) 

+ L 10 (U + F*UF L ^) + H 1 {F^F R,y,v + F L ^F L ^). (4.31) 

The meaning and renormalization of low-energy constants in Eq. (4.31) are independent of 
the symmetry breaking sector and, in particular, of the infrared dimension of the quark 
mass. The remaining i? -independent terms in Eq. (4.30), cf. £ 22 and C 04: , are absent from 
the expression for C^': with quark mass of dimension 2, these terms would count as 0(p 6 ) 
and 0(p 8 ) respectively. On the other hand, all terms but £40 contained in C^' are already 
included either in or in C {3) . Consequently, £ (2) + £ (3) + £ (4) not only encompasses all 
terms of the standard but, in addition, it contains new terms of the type £ ( - 3 - ) , £22 

and £o4- This phenomenon is general. Order by order, the improved contains the 
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standard perturbation theory as a special case: It contains more parameters and it could 
well fit the experimental data even when the standard fails. 

The one-loop contribution to the tt — n amplitude A(s\tu) has been worked out within 
the standard chiral perturbation theory in Refs. The result can be expressed in terms 

of four constants: a = (3 (close to 1), the shift <5a 4 (8a 3 = in this case) introduced in 
Eqs. (4.21), and two linear combinations of the renormalized constants Li,L2 and L3. In 
the improved x?T, the one- loop 0(p 4 ) amplitude contains, in addition, two parameters 
which arise from the new terms £22 and £04 in £ ■ Working with the amplitudes T, U, V 
(the formula (3.2) is valid up to and including two loops), one may obtain a closed form for 
the one-loop amplitude which encompasses both alternatives of chiral perturbation theory. 

Let <f^(s) denote the effective dimension-c? contribution to the real part of the partial 
wave amplitude f a (s), (a = 0, 1, 2), introduced in section III B: 

ite/a(a)=X>?>. (4-32) 

d>2 

From Eqs. (4.27) one finds 

<p { *\s) = 6P(s + K ) 

<p¥\s) = -3P(s + K 2 ), (4.33) 

where 

* S ^(-%-\)Ml (4.34) 

Similarly, the real parts at the 0(p 3 ) level are 

^ 3) = 5M 2 n 5a 3 , = 0, = -2K &*s, (4.35) 

where 5a 3 is given by Eq. (4.23). For d > 3, the real parts are no longer defined by the 
tree amplitude alone. The 0(p d ) contribution to the imaginary part of the partial wave 
amplitudes Im (s) can be expressed for s > AM 2 through elastic unitarity: 
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Jro/M( a ) = «/£_^£ £ ^\s)^\s). (4.36) 

» S d!+d 2 =d 

This result is an exact property of amplitudes for 4 < d < 8. 

The one-loop level contains d = 4, d = 5 and d = 6 contributions to the scattering 
amplitude A(s\tu). In the following, we shall merely concentrate on the leading 0(p 4 ) part. 
The corresponding components of the functions T, U, Vf\ will be denoted as t/^^s), U^l d (s), 
and V^ d (s). The discontinuities of these functions are given by the 0(p 4 ) absorptive parts 
Imf^\ following Eqs. (3.4). Hence, the 0{p 4 ) one-loop amplitudes T,U,V can be written 
as 

= l{[^o\s)} 2 + 2[^\s)f}L(s, /i 2 ) + a 4 (/i 2 ) + a^ 2 )s 2 
= ^{2[^ 2) (^)] 2 - 5[^ 2) ( S )] 2 }L( S ,/i 2 ) 

where L(s, jj 2 ) is the loop integral subtracted at the point s = — fj 2 : 



t _ s + ji 2 f°° dx 1 jx - 4M 2 



7T 74M2 X + /i^ X — S V X 

The constants a n (fi 2 ) and /3 n (/z 2 ) behave in the chiral limit as M™. They describe the most 
general polynomial part of T, U and which is 0{p A ) and takes into account the freedom 
(3.5). These constants represent renormalized tree contributions of the d — 4 part of £ e //. 
Their dependence on the subtraction point /i 2 can be determined by demanding that the 
scattering amplitude A(s\tu) be //-independent. Following Appendix B, this requirement 
is equivalent to the conditions 

9 r T«(s) = 5T(s), -^U^(s) = 5U(s), J^V^(s) = 5V(s), (4.39) 



dfj, 2 ' d\i 2 ' dfi 2 

where ST, 5U and 5V are of the general form (3.5). Taking into account the s-independence 

d 

of 7— — L(s,/i 2 ) and L(0, /z 2 ) = — L(— /U 2 ,0), the solution of Eqs. (4.39) can be easily found: 

6\l 2 



4 Notice that 0(p N ) terms in V contribute to the scattering amplitude A of Eq. (3.2) as 0(p N+2 ) 
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a (fi 2 ) = a (0) + 18p 2 L(-fi 2 ) 
P (fi 2 ) = /3 (0) 

a 4 (/i 2 ) = a 4 (0) + (lid 2 - 32(3 2 )M^L(-fx 2 ) 
&(/i 2 ) = ft(0) + 18/3(5d - 2/3)M 2 L(V) 



(4.40) 



In these equations we have denoted 

1 



L(s) = L(s,pL 2 = 0) 



7T 



aln( 



CT+ 1' 



1- 



4M3 



(4.41) 



In the following we shall work at jj 2 = 0. 

The constants ao and [3q are related to the low-energy parameters Li,L 2 and L 3 which 
occur in the expression (4.31) for £ 4 q. One gets 



«o(0) 

A)(o) 



^^[l: + + -l 3 - -K/^ 2 )l 

4ttF 4 L 1 2 2 J 4 v ^ ;J 

1 



8ttF 4 



(If 2 — — L 3 



10247r 3 F n 4 



(4.42) 



where 



32vr 2 



' Ml 1 M| 
/i^ 8 // 2 



+ 
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(4.43) 



and p 2 denotes the renormalization scale introduced in Ref. ||. The renormalized constants 
L[, L r 2 are /independent, whereas L 3 and L r 2 — 2L\ are not. Furthermore, the combination 
U 2 — 2L\ should be suppressed by the Zweig rule or in the large N c limit. Notice that the 
constant (5q is independent both of [i 2 and of p 2 . The interpretation of the remaining two 
constants a 4 and /3 2 depends on the effective dimension of the quark mass. In the standard 
chiral perturbation theory, these constants can be expressed in terms of the shifts of the pion 
mass and decay constant, as calculated within SU(2) x SU(2) perturbation theory In 
the improved x?T, ct 4 (0) and /3 2 (0) are independent parameters which describe respective 
contributions of new terms £04 and £22 in the 0(p 4 ) effective Lagrangian £ (4) . The explicit 
relationship between a 4 (0),/? 2 (0) and the low-energy parameters of CS^ is of no direct use 
in the present paper and it will be given elsewhere. 
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Concluding this section, it is worth noting that the low-energy theorem of Sec. Ill 
considerably simplifies the calculation of two-loop contributions to A(s\tu): For d < 8, all 
0(p d ) terms can be obtained by a straightforward combination of Eqs. (3.2) and (3.4) with 
the unitarity condition (4.36). Up to and including two loops, the expansion of the 
n — n scattering amplitude can be viewed as an iteration of the Roy-type Eqs. (3.7a) and 
(3.7b). The corresponding polynomials P a (s) appearing at a given order 0(p d ) are then 
defined in terms of the renormalized low-energy constants of the Lagrangians or , 
according to the effective dimension of the quark mass being respectively 2 or 1. 



V. DETERMINATION OF PARAMETERS OF £ eff 
FROM vr - Ti SCATTERING DATA 

Suppose one has enough experimental information to perform the program formulated 
in Sec. Ill and to reconstruct the low-energy amplitude A(s\tu). Let us call the result of 
this reconstruction A exp (s\tu) and the corresponding T, U, V amplitudes given by Eqs. (3.4) 
T eX p, U exp and V exp respectively. We would like to compare the experimental amplitude A exp 
with the theoretical amplitude A th given in Sec. IV in a whole low-energy domain of the 
s-t-u plane including the unphysical region. Such a comparison should lead to a detailed fit, 
which in turn should provide a rather precise determination of low energy constants entering 
A t h. In particular, we would like to measure the parameter a and, in this way, let Nature 
tell us whether it prefers a quark mass of effective dimension 1 or 2. The theorem proved in 
Sec. Ill considerably simplifies the above task: Neglecting 0(p 8 ) contributions, the equation 

A exp {s\tu) - A th (s\tu) = 0, (5.1) 

which is supposed to hold in a crossing symmetric domain of the Mandelstam plane, is 
actually equivalent to a set of three single-variable equations, 

T exp (s) -T th (s) = 5T(s) 
U exp (s) - U th (s) = 5U(s) 
32 



V exp (s) - V th (s) = 6V(s) (5.2) 

valid in an interval of s. The functions 5T(s), 5U(s) and 5V(s) are the arbitrary and irrel- 
evant polynomials given by Eq. (3.5). In this section, we will analyze Eqs. (5.2). Hereafter 
we systematically set = 1. 
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A. One-loop precision 



The functions T exp , U exp and V exp are given by Eqs. (3.6). Up to and including (one-loop) 
order 0(p A ), the theoretical amplitude reads 



(5.3) 



(and likewise for U and V), where the tree and leading one-loop contributions are presented 
in Eqs. (4.27) and (4.37) respectively. We shall concentrate on real parts of Eqs. (5.2). 
Let us denote the partial wave integrals appearing in Eqs. (3.6) as 

s 3 [ a2 dx Im f a (x) 



Ms) = -f 

7T J 



A X° X — S 



a = 0,2 



, . . s 2 r A2 dx 1 Im fi (x) 

Ms) = -+ 



(5.4) 



71 J A X 2 x — 4 x — s 

It is convenient to take linear combinations of the Eqs. (5.2) for T and U and isolate the 
contributions of / = and I = 2 s waves. The resulting equations can be written as 
(a = 0,2) 



(3 2 N % 
6tt 



S + K a ) 2 D(s) +p a (s), 



(5.5a) 



where 



iV = 36, 



No = 9 



(5.6) 



and (w 



4 

1 - - 

s 



1/2 



D(s) = 6tt ReL(s), 
D(s) = 12 + 6wln 



w 



s<0, s>4, 

1 + w 

D(s) = 12 - ^warctanw" 1 , < s < 4. 



(5.7) 



The p a (s) are two third-order polynomials, whose coefficients are given in terms of (i) three 
constants U [cf. the first of Eqs. (3.6)], (ii) the parameters a,(3,ao(0) and a 4 (0) defined 
in terms of £ e ff, an d (iii) the irrelevant five constants that characterize the polynomial 
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ambiguity (3.5). The explicit expression for the coefficients of p a (s) can be easily read off 
from Eqs. (3.6), (4.27) and (4.37). Similarly, the U-equation (5.2) can be written as 

Ms) = ^(s-A)D(s) + q(s), (5.5b) 

where q(s) is now a second order polynomial with coefficients given by linear combinations 
of three parameters Vi [cf. the last of Eqs. (3.4)], the £ e ff parameters (3 (0) and /3 2 (0) and 
the irrelevant constants yi. For small s, the function D(s) behaves as 

D(s) = s + ^ + Ois 3 ). (5.8) 

On the other hand, the functions a (s) and 0i(s) defined in (5.4) behave as 0(s 3 ) and 0(s 2 ) 
respectively. The polynomials p a (s) and q(s) should be such to insure this small s behavior 
on the right hand sides of Eqs. (5.5a) and (5.5b). Using (5.8), one easily finds 

Pa^) = ^ - ^ a (n a + 20)s 2 + r a s 3 | a = 0, 2 

^) = ^{-s( S -4) + r lS 2 }, (5.9) 

where to,ti,t 2 are three yet undetermined parameters. Eqs. (5.5a) and (5.5b) now take 
the form (a = 0, 2) 

B 2 N ( 1 1 

Ms) = ^ [s 2 D(s) + 2s[D(s) - s}n a + [D(s) - s - -s 2 ]K 2 a + r aS 3 | 

Ms) = ^{(s~ 4)[£>(a) ~s]+ r lS 2 } . (5.10) 

Once the experimental phase shifts are known, one can compute the integrals <j)(s) on left- 
hand side of Eq. (5.10) and fit them with the corresponding right-hand side. The parameters 
of the fit are a, (3, r , Ti, r 2 . At this stage, one does not need to know the subtraction 
constants tj and Vi in the dispersion relations (3.6). The latter are needed, however, if one 
wants to measure the four parameters of namely a (0), (3 (0),56t + a 4 (0) and /^(O). 
(Remember that the parameters c*o(0) and Po(0) determine the two linear combinations 
(4.42) of the low-energy constants L ± , L 2 and L 3 that appear in the £ 40 -part (4.31) of £ e ff-) 
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Indeed, comparing coefficients of polynomials on both sides of Eqs. (5.9), one gets 11 linear 
relations among the "experimental" constants t , t 2 , t 3 , vi, v 2 , v 3 , the four parameters of C e ff 
mentioned above, and the irrelevant five constants x,yo,yi,y 2 ,y3- Eliminating the latter, 
one can express the four C e ff parameters as 

P 2 

«o(0) = t 2 - T7- {2k (k + 20) + k 2 (k 2 + 20)} 

1U7T 

93 2 

A)(0) = ^ 2 + — (l-8r + 5r 2 ) 

7T 

33 2 

+ TFT {8«o(«o + 20) - 5k 2 (k 2 + 20)} , (5.11a) 

407T 

and 

5a + 0:4(0) = t — a — - — (2kq + k%) 

on 

P 2 (0)=v 1 -9P + ^ {54-84} 
63 2 

+ — {5(/€ 2 -2t 2 )-8(k -2t )}. (5.11b) 

7T 

The remaining two equations do not involve any parameter of C e ff to be determined. They 
read 

2 

t 3 = - f —{2r + T 2 ) 

7T 

9B 2 

V3 = ^(l-8To + 5r 2 -r 1 ). (5.12) 

The two Eqs. (5.12) should be merely expected to measure the strength of neglected 
two-loop and contributions, rather than represent a true constraint on the fit based on 
Eqs. (5.10). 



B. Fits to Roy-type equations (3.7a) and (3.7b) 

In order to reconstruct the amplitude A exp (s\tu), one needs a complete set of pion-pion 
phase shifts 5 a (s), (a = 0, 1, 2). (By complete we mean that they extend in energy from the 
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threshold to A < 1 GeV for all three isospins and are dense enough in the interval to allow 
adequate numerical evaluation of our dispersion integrals.) 

There exists only one complete set of pion-pion scattering phase shifts (extrapolated from 
experimental dataQ) that has been published in numerical form, namely that appearing in 
the paper of Froggatt and Petersen |T7|] . They provide values for S a (s) — without quoted 



errors — at 20 MeV energy intervals in 4M 2 < s < A 2 , for a = 0, 1, 2. The phase shifts S a 
come from an analysis following that of Basdevant et al. |||], which employs a truncated set 
of twice-subtracted Roy equations, makes a particular choice of parametrization for f a (fixing 
the 1 = scattering length, ao) and uses a Regge type model for estimating the high energy 
contributions to the dispersion integrals. Data were taken from the Estabrooks-Martin 
analysis |27] of the CERN-Munich experiment on tcN — > tttiN |28]. Although Basdevant et 



al. J23] present graphical results for several choices of values of Qq in their work, numerical 



results are only presented in the subsequent paper of Froggatt and Petersen p7|| , and only 
for the unique choice a = 0.3. 

We first check to what extent the Froggatt-Petersen phases satisfy the version of the Roy 
equations set forth in Section III B. To this end, we compute the integrals on the right hand 
side of Eqs. (3.7a,b), using the 5 a from Froggatt and Petersen. Calling the result Ref^ HS (s), 
we then determine the parameters ti,Vi by minimizing 

E URef a LHS (s.) - Ref a RHS (s*)?, (5-13) 

a i 

where Ref^ HS is the real part of f a determined directly (via unitarity) from 8 a . This is not 
a proper % 2 fit, since no uncertainties can be included; consequently, no uncertainties can 
be quoted for the resulting constants. We find, however, that the values for experimentally 
determined constants are stable for reasonable variations in the energy interval used for the 
fit (see Table [I]). The fit over the largest range, 4 < s < 25, is excellent: Ref^ HS and 
Refa HS agree to 1% over nearly all the interval, the sum in Eq. (5.13) being O(10~ 4 ) for 63 



'For a recent review of experimental ir — tt scattering data, see [25,2' 
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data points. We see no need to present the results graphically: Ref^ HS and Ref^ HS would 
be indistinguishable. Instead, the values of Ref^ HS and Ref^ HS are compared in Table f|, 
for 21 energies included in the sum (5.13). We thus conclude that the Froggatt-Petersen 
phases indeed give a solution of our set of triply-subtracted Roy equations, for the values of 
parameters ij and Vi summarized in Table |. (Notice that the parameters £3 and v 3 are poorly 
determined, but that they are sufficiently small not to affect the analysis at the 0(p 4 ) level.) 
The corresponding low-energy amplitude A exp (s\tu) will be confronted with the theoretical 
prediction A th shortly. 



-K" e 4-decay experiments |29j are consistent with the value ao = 0.30 for the scattering 
length, characteristic of Froggatt-Petersen phases, but standard predicts a lower value, 
namely ao = 0.20 ± 0.01 ||. It would be desirable to have complete sets of phase shifts that 
fit both experiment and Roy equations for other values of ao < 0.30. These are not available.^ 
For this reason, we must use ad hoc extrapolations down to threshold of existing data at 
energies E > (500 — 600) MeV obtained from tcN — > tcttN and tcN —>■ 7T7tA production 
experiments. One such extrapolation has been recently considered by Schenk 0] using a 
simple parametrization 

7 ,s - 4 X ,5-4, 



tan. 5 As) 



a t + b l (^) + c l (^) 2 
4 



b i = b i -a i — — + (aif. (5.14) 
s - 4 

for the two s waves (i = 0.2) and a similar formula for the p wave. The scattering lengths 
^ and the slope parameters 6j are fixed at their values predicted by the standard one-loop 
XPT §: 

a = a° Q = 0.20, a 2 = a 2 = -0.042, a x = a\ = 0.037 

b = b° = 0.24, b 2 = b 2 = -0.075 . (5.15) 

The remaining parameters are determined by fitting the data obtained from various analyses 



J.L. Basdevant, private communication. 
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of dipion production experiments ||28|| . For the I = s wave, Schenk uses the Ochs energy- 
independent analysis]] of the CERN-Munich experiment [p7|] , covering the energy range 610- 
910 MeV. For his best fit - called solution B - no x 2 or error bars are quoted. Instead, two 
additional sets of parameters Co and sq = Eq called "A" and "C" are given that bracket 
together both the Ochs data and the well-known data by Estabrooks and Martin \27\. A 
similar procedure is adopted for the I = 2 s wave, whereas the parameters of the p wave 
are determined from the experimental p mass and width. Results of this analysis and more 
details can be found in Ref. fl3~0~| . 

In this way, the parametrization (5.14) provides a complete set of phases - hereafter 
referred to as Schenk B - that fits the data at higher energies and uses the threshold 
parameters (5.15) of the standard ;\RT. Using this set, we have performed exactly the same 
kind of fit to the Roy-type Eqs. (3.7a) and (3.7b) as in the case of Froggatt-Petersen phases. 
Surprisingly enough, we find this fit at least as good as in the case of the Froggatt-Petersen 
phases, despite the fact that the Schenk B phases were not obtained using Roy equations 
or any other crossing-symmetry correlation among the three lowest partial waves.Q The 
resulting parameters ti and V{ are given in the second half of Table |, and the quality of 
the fit can be appreciated from Table [TV]. Unfortunately, we do not see any simple way to 
associate the Schenk B phases, and the corresponding parameters ti and Vi, with a set of 
errors which would be deduced from statistical errors of the experimental data used at the 
beginning and which would respect the correlations imposed by the Roy equations. The 
same remark applies to the set of phases of Froggatt and Petersen. 



7 The data by Ochs can be found in his unpublished thesis [31]. We are indebted to Dr. J. Gasser 
for communicating these unpublished data to us. For the reader's convenience, they are reproduced 



in our Table III. 



? A. Schenk, private communication. 
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C. Determination of parameters a, (3, L\, L<i and L3 
from a complete set of phase shifts 



The next step is to confront the empirical amplitude A exp with the amplitude A t h com- 
puted from chiral perturbation theory. In particular, the two solutions of the Roy-type 
equations (3.7a) and (3.7b) described above can be used to measure the parameters a, (3 
and, through Eqs. (4.42), two linear combinations of the low-energy constants L\,U 2 and 
L3, defining the four-derivative terms in C e ff- The measurement is based on Eqs. (5.10). 
First, one evaluates the three functions 4> a (s), a = 0,1,2, defined in Eqs. (5.4), using the 
complete sets of phase shifts exhibited in Section V B. The results are represented graph- 
ically by continuous lines in Figs. 3a,b,c for the case of Froggatt-Petersen phases and in 
Figs. 4a,b,c for the Schenk B set. Next, one fits the "experimental" functions <p a (s) with 
the theoretical expression represented on the right-hand side of Eqs. (5.10). The parameters 
of the fit are a, f3 and r , T\, t 2 . [Recall that the K a are defined in terms of the ratio ot/ (3 - 
see Eq. (4.34).] 

The range in s in which the fit is performed should not exceed the range in which the 
0(p 4 )-order %PT may actually be expected to apply. On the other hand, this range should 
be large enough to permit a sensitive determination of parameters. For this reason, it might 
be misleading to consider exclusively the physical region s > 4 Pl,|2"6]|. In the following, we 
use the interval —4 < s < 8, which most likely represents a rather conservative choice. 

From Figs. 3a,b,c and 4a,b,c one observes a large difference in scale of individual a : 0o 
is typically an order of magnitude or more larger than 02 and nearly two orders of magnitude 
larger than 0i. For this reason, we first fit the function O , determining the three parameters 
a, (3 and To- Then, using the values of a and j3 obtained in this way, we perform two single- 
parameter fits to 02 and 0i, determining T2 and T\ respectively. In the absence of error bars 
for (f> a {s), it is impossible to perform a true \ 2 fit- Instead, we minimize the sum of squares 
of the difference between the left- and right-hand sides of Eqs. (5.10), for 66 equidistant 
points in the interval —4 < s < 8, giving the same weight to each point. 
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In all cases, the parameter (3 = F%/F 2 should remain close to 1, and the fit should be 
constrained by this condition. We require 

/3 < 1.17, (5.16) 

corresponding to the lower bound F > 86 MeV. This bound is consistent both with existing 
standard ;\PT estimates and with the improved xPT formula (4.20). Leaving the ratio 
a/ (3 unconstrained in the minimization procedure, one tests - for a given set of data - the 
relevance of the improved %PT. The corresponding fits are represented by dashed curves in 
Figs. 3a,b,c and 4a,b,c. The corresponding best values of the parameters are 

a//? = 4.20, ,3 = 1.17 
r = -0.263, n = 3.75, r 2 = -0.540 (5.17a) 

for the set of Froggatt-Petersen phases, and 

a//3=1.63, (3=1.17 
r = -0.032, 71 = 3.68, r 2 = -0.640 (5.17b) 

for phases of the Schenk B set. On the other hand, in order to test the compatibility of the 
0(p 4 ) standard ;^PT with a given set of data, one further restricts the fit by requiring 

a = /3<1.17. (5.18) 

Results of the minimization with this constraint are represented by dot-dashed curves in 
Figs. 3 and 4. The best values of parameters corresponding to this constrained fit are 

a = (3= 1.17, t = -0.414, n = 3.75, r 2 = -0.661 (5.19a) 

and 

a = (3 = 1.17, r = -0.045, n = 3.68, r 2 = -0.653 (5.19b) 

for the Froggatt-Petersen and Schenk B sets of phases respectively. 
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A few remarks are in order. The Froggatt-Petersen data are considerably better fit 
in terms of a larger value (5.17a) of the ratio a/ (3 than the standard %PT would permit, 
although without a true x 2 At we cannot be quantitative about this observation. The failure 
of standard %PT to describe the Froggatt-Petersen s wave is also apparent in Fig. 3a 
(dot-dashed curve). Concerning the p wave, the fit is reasonably good for both cases (Fig. 
3b), reflecting the fact that the theoretical calculation of 0i(s) senses the effective infrared 
dimension of the quark mass starting only at the two-loop level. It is worth noting that 
the best value a//3 = 4.20 is overcritical by 5%. This means that the Froggatt-Petersen 
1=0 s wave would be compatible with the vanishing of the qq condensate -Bo-0 From this 
point of view, the set of Froggatt-Petersen phases with = 0.30 appears as an extreme 
alternative. The opposite extreme is represented by the Schenk B set of phases. Since the 
latter incorporates a priori the values of scattering lengths and effective ranges as predicted 
by the standard %PT, it is not surprising that the corresponding best value for a/ (3 (5.17b) 
is considerably closer to 1 than in the Froggatt-Petersen case. Furthermore, Fig. 4a seems 
to indicate that, although the best value for a/ (3 is still as large as 1.63, this fact need 
not be significant. In the absence of error analysis, it is hard to be too affirmative in the 
interpretation of the Schenk B fit. 

It remains to exploit the additional information (values of constants t and v as well 
as the constants r resulting from our fits), in order to measure certain parameters of the 
dimension-4 component of £ e ff- Here, we merely concentrate on the constants L 1; L 2 and L 3 
characteristic of £40, Eq. (4.31), whose meaning and renormalization do not depend on the 
effective dimension of the quark mass. For this purpose, we have to determine the constants 
ao(0) and /3q(0) given by Eqs. (5.11a). Using the central values of the parameters t 2 and v 2 
(the second column of Table |) and the best values for a/ (3 and the r's, as determined in 
the previous fits, one gets a (0) = 5.81 x 10~ 4 , (3q(0) = 3.99 x 10 -3 for the Froggatt-Petersen 



3 This critical case has been considered earlier [32] 
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solution, and a (0) = 5.87 x 1(T 4 , #,(0) = 2.07 x 10" 3 for the case of Schenk B phases. These 
numbers are easily converted into information on the constants £1,2,3; using Eqs. (4.42) and 
(4.43). Assuming the Zweig-rule (or large- N c ) relation L T 2 — 2L\ = 0, and identifying the 
running scale Jx with the 77 mass, as done in Refs. P,|33|, one obtains 



L\ = 2L\ = 1.34 x KT 3 , L 3 = -4.50 x 10~ 3 (5.20a) 

for the Froggatt-Petersen data, and 

L r 2 = 2L\ = 0.56 x 10~ 3 , L 3 = -2.15 x 10~ 3 (5.20b) 

for the set of Schenk B phases. It is gratifying to see that these values - especially (5.20a) - 
compare well with other determinations based on standard xPT P,|33|. Indeed, there is no 
reason why the purely derivative terms in C e ff should be affected by questions concerning 
the symmetry breaking sector. 

D. Estimates of errors in the direct measurements of a/ (3 

The uncertainties in the values of the parameters a,/3 and r a arise from uncertainties 
in the functions cb a ; these uncertainties, in turn, arise from uncertainties in the phase shifts 
S a over the range of integration in Eqs. (5.4). As we have noted, there is no set of phase 
shifts S a which exists, together with corresponding errors, in this energy range. In the 



present subsection, we extend the extrapolation method of Schenk ^0] , described above, to 
construct several sets of 1=0 phase shifts Sq, together with estimated errors, in the necessary 
energy interval. In this way, we obtain values and estimated errors for the parameters a/ (3 
and To for each extrapolated data set. Only 1 = phases are considered. In fact, we could 
treat 1=1 phases similarly (although the insensitivity of <pi to a makes this relatively 
uninteresting); in any case, the paucity of experimental data on / = 2 makes the production 
of a complete set of phase shifts impossible without a more extensive recourse to the use of 
the Roy equations, as in the analysis of Basdevant et al. [231 . 
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The two original sets of phase shifts (with corresponding errors) used are that of Ochs and 
that of Estabrooks and Martin. These were each obtained independently from analysis of 
the same CERN-Munich experiment. The first step is to extrapolate Sq down to threshold, 
using the Schenk formula (5.14). The Ochs phases are fit over the energy range 610-910 
MeV, i.e., using all his data for which no inelasticity is suggested. (See Table |TJ.) The 
Estabrooks-Martin phases are fit over the energy range 570-910 MeV, i.e., using all their 
points in the elastic scattering region except for their first three lowest-energy points, which 
appear to be less trustworthy. In performing the extrapolation, the scattering length do is 
fixed and the remaining parameters bo, cq, Eq are determined by minimization of x 2 using 
the phases and errors given to us. We show in Fig. |5] the results of this fitting procedure 
for the choices a = 0.20 (preferred by standard %PT) and 0.26 (preferred by if e4 -decay 
experiment) for each of the two data sets; the resulting parameters are given in columns 
2-5 of Table |V|. The x 2 f° r these fits is quite good. (We note in passing that the data of 
Estabrooks and Martin is not well described by the parameters ao = 0.20, b = 0.24 which 
characterize the Schenk B solution.) The next step is to estimate the uncertainty in the 
extrapolated phases Sq. Since the dominant parameter (after a , which is fixed) is bo, and in 
view of the strong correlations among the parameters, we proceed as follows: for fixed values 
of bo larger than its value for Xmim the minimum-x 2 value, fit the data by allowing cq and 
E Q to vary freely, and find the values of b , c , E which give x 2 — Xmin + lj ca ^ this solution 
"a", in analogy with Schenk's notation; repeat this procedure for fixed values of b smaller 
than that for Xmini can this solution "c"; the uncertainty in the phase shift S (Ei), for each 
value of E = E{, is then estimated by interpreting the variation of Sq{E{) from its solution 
"a" value to its solution "c" value as ±1 standard deviation in 5q{Ei). (This is similar to 
the procedure adopted by Schenk, although he allows much greater variation - leading to 
much larger uncertainties - in order to bracket both Ochs and Estabrooks-Martin phases at 
the same time.) 

Now, for each of the four sets of phase shifts So, obtained by the extrapolation procedure 
described above, we may make the comparison for O as done for the Froggatt-Petersen and 
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Schenk B phases in Section IV C. However, we now have the important advantage that a 
true x 2 fit is possible, so we can have some idea of the precision with which the resulting 
parameters are determined. For each set, we make two fits: one, corresponding to standard 
xPT, for which we fix a = f3 < 1.17; the other, corresponding to improved xPT, for which 
(3 < 1.17 but a is allowed to vary freely. The fits are all performed over the same interval 
— 7 < s < 9. Results of the determination of the parameters a/ (3 and To are given in columns 
6-9 of Table [V|; the reader can judge the quality of the fits from the plots of 0o given in 
Figs.^a,b,c,d. The solid curves represent the parametrization of improved xPT, while the 
dashed curves represent that of standard xPT. It is clear, both from the large x 2 values 
tabulated for the standard xPT fits and from examination of the dashed curves in Figs.|6| 
that standard xPT is not compatible with these phase shifts. For this reason, we quote no 
result for r for this case. On the other hand, improved %PT can easily accommodate such 
data. It is important to note that, for a given set of phase shifts, the parameters a/ (3 and 
t are very well determined by the improved %PT fit. 

As a check on our procedure of estimating errors, we have also used a more "conservative" 
procedure, viz., vary all non-fixed parameters within their one-standard-deviation limits 
to produce solution "a cons " , taking max(6 ),max(co),min(s ), and solution "c cons ", taking 
min(6o),niin(co),max(so); then compute the "conservative" uncertainties in the phase shifts 
S (Ei) using these solutions as we did for solutions "a" and "c" before. Clearly, this method 
does not take into account the strong correlations in b , c , s . Thus, when the consequent 
phase shift errors are used in the fitting of <po, these larger errors result in larger errors in the 
experimental function O - The result is then a x 2 roughly half of that previously obtained, 
and errors in a/ (3 and To roughly 2-5 times larger. Nevertheless, the best fit is the same. 

VI. SUMMARY AND CONCLUSIONS 

A new framework for testing the convergence rate of chiral perturbation theory is pro- 
posed. One first replaces the standard expansion of the effective Lagrangian by a more 
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general expansion that is as systematic and unambiguous as the standard %PT. In addition 
to the usual terms, the new expansion involves at each given order new contributions that 
the standard ;\RT relegates to higher orders. The size of these additional contributions can 
then be tested experimentally, in particular in low-energy n — n scattering. Unless these 
contributions turn out to be small, the improved xPT has, in principle, more chance to 
produce a rapidly convergent expansion scheme. 

A new low-energy theorem is presented which provides the general solution of constraints 
imposed by analyticity, crossing symmetry and unitarity on the tt — n scattering amplitude, 
neglecting 0(p 8 ) contributions. Applications of this theorem are threefold: 

i) First, it considerably simplifies the evaluation of the perturbative 7r — tt amplitude up 
to and including two loops. This applies both within the "standard xPT" and within the 
more general "improved xPT" which contains the former as a special case. In both cases, 
the calculation reduces to the iterative insertion of the unitarity condition (4.32) into the 
dispersive integral for the functions T,U and V in Eq.(3.2). The improved xPT one-loop 
amplitude is worked out in detail in Sec. IV. The two-loop amplitude can be easily calculated 
along the same lines. The reason why the formula (3.2) no longer holds beyond two loops 
resides in new 0(p 8 ) effects in the absorptive part: inelasticities and higher partial waves. 

ii) Next, the low-energy theorem of Sec. Ill can be used to constrain the low-energy scat- 
tering data and to fully reconstruct the corresponding amplitude. The formula (3.2) implies 
a particular truncation of the infinite system of Roy equations, under a rigorous control of 
chiral power counting: Neglected contributions are 0(p 8 ), whereas in the original form of 
the Roy equations p(| the model-dependent "driving terms" are of the same order 0(p ) as 
the effects we are looking for. A complete set of low-energy phases $q,5q and <5}, together 
with the six subtraction constants t and v for which the Roy-type Eqs. (2.7a) and (2.7b) 
are satisfied to a reasonable accuracy (see Tables [IT] and [TVD, define up to 0(p 8 ) corrections 
the scattering amplitude A(s,t,u) in a whole low-energy region of the Mandelstam plane 
including the unphysical region. Two examples of such a complete low-energy amplitude 
are given, based on phase shifts published by Froggatt and Petersen |P7| and by Schenk |3(| 
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respectively. They are both compatible with existing irN — * ttttN and K e4 experimental 
data. 

iii) Finally, the low-energy representation (3.2) simplifies the direct comparison of the 
perturbative amplitude A th (s,t,u) with the amplitude A exp (s, t, u) reconstructed from the 
data. In particular, parameters of C e ff contained in A th can be measured through a detailed 
fit of the amplitude A exp (s, t, u) over a sufficiently large portion of the Mandelstam plane in 
which the low-energy expansion can still be taken as valid. The fit is particularly sensitive to 
the ratio a/ (3 which parametrizes the leading 0(p 2 ) amplitude. The improved requires 
1 < a/ (3 < 4, whereas the special case of the standard xPT corresponds to a/j3 = 1. The 
ratio a//3 is related to the value of the QCD parameter 2mB in the units of pion mass 
squared and, via the pseudoscalar mass spectrum, to the quark mass ratio r = m s /m. 

Examples of measurement of a//3 exhibited in this paper illustrate the lack of sufficiently 
precise experimental information on low-energy n — n scattering. For a fixed value of the 
scattering length Oq, the statistical errors of the production data on 5° are estimated to show 
up as errors in the measured values of a/ (3 of the order of a few percent. On the other hand, 
for different values of a° in the experimental range a$ = 0.26 ± 0.05, and for different sets of 
production data, the resulting values of a/ft vary between 1.5 and 4.2. The two complete 
low-energy amplitudes mentioned above correspond to these two extremes. In particular, 
the Froggatt-Petersen phases (for which Oq = 0.30) are compatible with the vanishing of the 
condensate B and with the critical value of the quark mass ratio r = r\ ~ 6.3. 

The suspicion that a bad convergence of the standard might bias the usual con- 
clusions that r = m s /m ~ 25.9 and 2mB ~ M% is at least well motivated, but clearly 
it requires confirmation. In order to produce a truly unbiased measurement of these fun- 
damental QCD parameters, the method developed in this paper can prove useful provided 
that it is supplied with more accurate experimental information on low-energy n — n phase 
shifts. The current imprecision, illustrated by error bars as large as in a° = 0.26 ± 0.05 can 
hide all cases of interest, including the intriguing critical case (qq) = 0. Here, one faces a 
challenge of fundamental high precision low-energy experimental physics. 
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APPENDIX A: NOTATION AND CONVENTIONS 



In the first part of this appendix, we fix the notation and normalization for the scattering 
amplitude. We then exhibit the main properties of the crossing matrices C [Eq. (3.9)]. 

The S- matrix element for the transition 7r a + 7r 6 — > n c + -K d , where a, 6, c, d are pion isospin 
indices, is connected to the T-matrix element by the relation: 



(cd\S\ab) = (cd\ab) + i(2n) 5 (p a + p b - p c - p d ) T abM 



(Al) 



The T-matrix element can be written in terms of isospin invariant amplitudes; taking cross- 
ing symmetry into account, the decomposition reads: 



T a b,cd{s,t,u) = A(s\tu)5 ab 5cd + A(t\su)5 ac S bd + A(u\ts)5 ad 5 bc , 
where s, t and u are the Mandelstam variables, 

S = (Pa + Pbf, t = (p a - Pc) 2 , U={p a - p d ) 2 . 



(A2) 



(A3) 



The amplitude A(s\tu) is symmetric in the variables t, u. (The amplitude A(t\su) is obtained 
from A(s\tu) by the exchange of variables s, t and by subsequent analytic continuation.) 
The s-channel isospin amplitudes [Eq. (3.8)] are related to the amplitude A by 

(z i i w N 

1 



I F (o)\ 



F W 



(s,t,u) 



32tt 



3 1 
1 
1 



A(s\tu) 
A(t\su) 
A(u\ts) 



(A4) 
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The partial wave expansion is 



= £(2€+l)P,(cos0)//(a), 



(A5) 



where s = 4(M^ + g 2 ),t = —2q 2 (l — cos9). With this normalization, the elastic unitarity 
condition for // takes the form: 



Imfl{s) 
//(«) 



s-AMl 



JSUs) 



sin 5p(s). 



The crossing matrices C [Eq. (3.9)] have the following forms: 



and satisfy the relations 



a 



tu 



n 2 — n 2 



1, 



(A6a) 
(A6b) 





/ 1/3 


1 


5/3 ^ 









o] 




' 1/3 


- 1 


5/3 ^ 


C s t — 


1/3 


1/2 


-5/6 


) C*u — 





- 1 





c 

5 ^ su 


-1/3 


1/2 


5/6 




V V3 


-1/2 


1/6, 













v 1/3 


1/2 


1/6, 



(A7) 



(A8) 



/^y fi _ fi fi fi fi 

^st^su — ^tu^st — ^su^tui 

ft fi _ fi fi fi fi 



(A9) 



It is worthwhile to notice that Eqs. (A9) imply that the eigenvectors A± with eigenvalues 
±1, respectively, of the matrix C tu satisfy 



C SU C s tA± — ±Cst,A 



(A10) 



These relations are extensively used throughout the calculation of Sec. Ill C. 

The invariance of the amplitude A(s\tu) under the exchange of the variables t, u permits 
us to construct rather easily all independent crossing symmetric polynomials of a given 
degree in s,t,u. It is convenient to take the variables t,u as independent. Since there are 
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k n = [n/2] + 1 



(AH) 



symmetric monomials in t, u that are homogeneous of degree n, the number of independent 
parameters in a general crossing symmetric polynomial of degree N is 

N 

K N = J2 kn- (A12) 

n=0 

Hence, the most general polynomial of degree three contains six parameters, as claimed in 



APPENDIX B: AMBIGUITIES OF THE AMPLITUDES T, U, AND V 

In this appendix, we determine the general expression for the transformations T — > 
T + 5T, U — > U + 5U, and V — > V + 5V that leave invariant the scattering amplitude A 
[Eq. (3.2)]. It follows that the variations 5T, 5U and SV must satisfy the equation 



In order to solve Eq. (Bl), one notices that only two out of the three variables s,t,u are 
independent. By successive differentiation with respect to independent variables, one obtains 
a set of simpler equations which can be solved easily. To simplify notation, let us define 



We first consider s and t as independent variables, and differentiate Eq. (Bl) first with 
respect to s and then with respect to t. We thus obtain the following two equations, where 
the primes indicate differentiation with respect to the arguments of the functions: 



Sec. Ill C. 




(Bl) 



f = 5T, g = 5U, h = 5V. 



(B2) 





(B4) 
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We now consider t and u as independent variables and differentiate Eq. (B4) with respect 
to t, obtaining the result 



h"(t) - h"(u) = 0, (B5) 

which indicates that h" is a constant, and therefore h is a quadratic polynomial, 

h{t) = ]-at 2 + bt + c, (B6) 

where a, b, c are constants. Using this result for h in Eq. (B4), we find a relation between / 
and g: 

f(u) = l -g{u) + \u 3 - 1(6 + 2Mla)u 2 + du + e, (B7) 

where d and e are constants. We then return to Eq. (B3), consider s and u as independent 
variables, and differentiate with respect to s. This implies 

f"(s) + 2 -g"{s) - |& - jja(4A4 2 - a) = 0, (B8) 

which becomes, after replacing / in terms of g [Eq. (B7)], 

g"(s) + as- 2(6 + 2M*a) = 0, (B9) 

the solution of which is 

g(s) = ~s 3 + (b + 2M*a)s 2 + ks + £, (BIO) 

where k and £ are constants. The expression for / then becomes 

/(*) = (d+|)a + (e + |). (Bll) 

Finally, upon substituting the expressions for /i [Eq. (B6)], g [Eq. (BIO)] and / [Eq. (Bll)] 
in the original equation (Bl), we obtain two constraints on the constants 

f AM 2 k 

( e +f) = -nr( d +f). (Bi2) 

16 

c = -(H4M^ + -M T 4 a). (B13) 
o 
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After relabeling the constants as follows, 



y =£, yi = k, y 2 = (b + 2M 2 a), y 3 = -a/6, x = d + fc/3, (B14) 

the functions /, g and h, and hence 5T, 5U and 5V [Eqs. (B2)], take the forms given in Eqs. 
(3.5). 

APPENDIX C: POLYNOMIALS AND KERNELS OF THE ROY TYPE 

EQUATIONS 

In this Appendix we explicitly list the polynomials P a (s) and the kernels W a b which 
appear in the Roy-type dispersion relations (3.7). 

Po(s) = 5t + j^ 2 [3s 2 + 2(s - AM 2 ) 2 ] + ^t 3 [2s 3 - (s - AM 2 ) 3 ] 
+ 2 - Vl (3s - AM 2 ) ' ^v 2 (s - M 2 )(s - AM 2 ) 

+ ±v 3 (5s - AMl){s - AMlf (CI) 
P 2 (s) = 2t + 2 -t 2 [3s 2 + 2(8 - AM 2 ) 2 ] + |t 3 [2 S 3 - (s - AM 2 ) 3 ] 

- ± Vl (3s - AM 2 ) + ^v 2 (s - M 2 )(s - AM 2 ) 

- ±v 3 (5s - AMl){s - AMlf (C2) 
P l (s) = 1 -(s-AM 2 )(v l + v 2 s) 

+ ^(s - 4M 2 )v 3 [s 2 - l( a - 4M 2 )(11 S - AM 2 )] (C3) 



f A (IT 4 /" 

- / -E^( S ^)/m/ b (a;) = -( S -4M w 2 ) 2 ( S -2M2) / 

7T J4M2 X b=0 71 J A 



AM% x 3 (x - AMI) 



- — (a - 4M 2 ) 3 / ^{Jm / (x) + 57m / 2 (x) 

07T J4M2 X 4 



1 M 2 drr Tr . . „ . , 2. ,,,,,,, ^„,9s [ a * dx Im fAx) 
~ —J2 W ^,x)Imf b (x) = —(s-AM^) 2 (s~2M 2 ) 

7T J AMI X 7T J 4 



6=0 „ 7 4 M2 x 3 (x - AMI) 
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-— {s-AMlf / ™{2Imf (x) + Imf 2 (x) 
12n J ami x 

\ C f t ^ *>*• AW = > - 4MJ)'(. - 2M?) j£ 

+ 1 (. - ,Mlf f * P - U + t^Ml) G (i^)] 

127T J4M2 X 3 \ X J \ X J 

x{2Imf (x) - 5tof 2 (x) + 9 + to h(x)} (C6) 

In Eqs. (C4)-(C6), the function G is defined as follows: 

r 1 y 3 4 2 4 4 
C( ) Jo + 3x x 2 ~^ x 3 x 4 ^ ~^ ^' 

G(0) = 1. (C7) 
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FIGURES 

FIG. 1. The region of allowed values for the ratios a/ (5 and m s /rh lies between the two curves 
shown. 

FIG. 2. Phase shift 5q — 8\ from data sets of Froggatt and Petersen [17| (dashed curve) and 
Schenk B |3fJ (solid curve) compared with experimental data [^J] from ET e 4-decay. 

FIG. 3. The functions 4> a (shown as solid curves) for (a) a = 0, (b) a = 1 and (c) a = 2, 
using experimental phase shifts given by Froggatt and Petersen |i~7fl . Comparison is made with 
theoretical fits: those of the standard xPT are shown as dot-dashed curves, while the improved 
xPT fits are shown as dashed curves. 



FIG. 4. The functions 4> a (shown as solid curves), using phase shifts from the Schenk B [3C 



parametrization of the phase shifts of Ochs [31]. The meaning of the curves is the same as in Figure 
3. 

FIG. 5. (a) Schenk-type parametrization of phase shift data of Ochs, fixing ao = 0.20 (solid 
curve) and ao = 0.26 (dashed curve); (b) Schenk-type parametrization of phase shift data of 
Estabrooks and Martin, fixing ao = 0.20 (solid curve) and ao = 0.26 (dashed curve). Details are 
given in the text. 

FIG. 6. The function <j)Q (shown as points with error bars), using Schenk-type parametrization 
of phase shift data: (a) data of Ochs, fixing ao = 0.20; (b) data of Ochs, fixing ao = 0.26; (c) data 
of Estabrooks and Martin, fixing ao = 0.20; (d) data of Estabrooks and Martin, fixing ao = 0.26. 
In each case, the solid curve represents the parametrization of improved xPT, while the dashed 
curve represents that of standard xPT. 
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TABLES 





TABLE I. Parameters 


resulting from fitting Eqs.(3.7a,b) 




Parameter 




Energy Range for Fit (MeV) 






300-580 


300-640 


300-700 


Using phase shifts of Froggatt and Petersen: 


to 


0.0206 


0.0207 


0.0208 




6.4 xl0~ 4 


6.5 xl0~ 4 


6.7 xl0~ 4 


h 


5.2 xl0~ 6 


3.5 xl0~ 6 


1.6 xl0~ 6 


Vl 


0.0764 


0.0760 


0.0755 


V2 


0.0021 


0.0020 


0.0020 


V3 


-1.3xl0~ 5 


-4.4xl0~ 6 


+3.5xl0~ 6 


Using phase shifts of Schenk, solution B: 


to 


0.0067 


0.0065 


0.0063 




4.6 xl0~ 4 


4.8 xl0~ 4 


5.1 xlO" 4 


h 


1.4 xl0~ 5 


9.9 xl0~ 6 


6.9 xl0~ 6 


Vl 


0.0697 


0.0695 


0.0693 


V2 


0.0021 


0.0020 


0.0020 


V3 


-8.4xl0~ 7 


2.3xl0~ 6 


6.5xl0~ 6 
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TABLE II. Comparison of the left and right hand sides of Eqs.(3.7a,b), using phase shifts of 



Froggatt and Petersen 



1=0 


1=1 


1=2 


ENERGY 


LHS 


RHS 


LHS 


RHS 


LHS 


RHS 


300. 


0.342 


0.344 


0.005 


0.006 


-0.029 


-0.029 


320. 


0.377 


0.376 


0.014 


0.012 


-0.043 


-0.041 


340. 


0.414 


0.411 


0.022 


0.019 


-0.055 


-0.052 


360. 


0.447 


0.445 


0.028 


0.027 


-0.067 


-0.064 


380. 


0.479 


0.477 


0.039 


0.037 


-0.080 


-0.076 


400. 


0.509 


0.507 


0.049 


0.047 


-0.090 


-0.087 


420. 


0.534 


0.534 


0.058 


0.058 


-0.100 


-0.099 


440. 


0.557 


0.556 


0.072 


0.071 


-0.113 


-0.110 


460. 


0.573 


0.572 


0.088 


0.086 


-0.122 


-0.121 


480. 


0.584 


0.585 


0.105 


0.103 


-0.132 


-0.132 


500. 


0.590 


0.591 


0.123 


0.122 


-0.142 


-0.143 


520. 


0.590 


0.591 


0.146 


0.144 


-0.152 


-0.153 


540. 


0.585 


0.586 


0.171 


0.170 


-0.161 


-0.163 


560. 


0.574 


0.576 


0.199 


0.199 


-0.171 


-0.172 


580. 


0.558 


0.558 


0.234 


0.233 


-0.178 


-0.181 


600. 


0.537 


0.538 


0.272 


0.273 


-0.186 


-0.189 


620. 


0.512 


0.513 


0.318 


0.319 


-0.194 


-0.196 


640. 


0.483 


0.484 


0.371 


0.371 


-0.201 


-0.203 


660. 


0.450 


0.450 


0.428 


0.429 


-0.209 


-0.209 


680. 


0.414 


0.413 


0.486 


0.490 


-0.214 


-0.214 


700. 


0.375 


0.373 


0.532 


0.531 


-0.222 


-0.217 
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TABLE III. Data from energy-independent analysis of Ochs' 



Energy (MeV) 


<5o (degrees) 


610 


56.3 ± 3.2 


630 


59.5 ± 2.9 


650 


65.6 ± 3.2 


670 


62.5 ± 3.5 


690 


68.8 ± 3.6 


710 


74.5 ± 3.8 


730 


79.4 ± 3.6 


750 


81.2 ± 5.7 


770 


79.9 ± 3.9 


790 


77.5 ± 5.7 


810 


84.1 ± 3.3 


830 


84.4 ± 2.6 


850 


87.1 ± 2.5 


870 


89.2 ± 2.5 


890 


93.2 ± 2.9 


910 


103.3± 3.2 


a Ochs, Ref. [31 
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TABLE IV. Comparison of the left and right hand sides of Eqs.(3.7a,b), using phase shifts of 



Schenk, solution B 



1=0 


1=1 


1=2 


ENERGY 


LHS 


RHS 


LHS 


RHS 


LHS 


RHS 


300. 


0.236 


0.234 


0.006 


0.005 


-0.053 


-0.052 


320. 


0.274 


0.272 


0.012 


0.011 


-0.064 


-0.063 


340. 


0.314 


0.312 


0.019 


0.018 


-0.075 


-0.074 


360. 


0.356 


0.355 


0.027 


0.026 


-0.087 


-0.085 


380. 


0.398 


0.397 


0.035 


0.034 


-0.098 


-0.097 


400. 


0.439 


0.440 


0.045 


0.044 


-0.109 


-0.108 


420. 


0.479 


0.479 


0.056 


0.055 


-0.121 


-0.119 


440. 


0.515 


0.515 


0.068 


0.067 


-0.132 


-0.130 


460. 


0.545 


0.546 


0.083 


0.081 


-0.143 


-0.142 


480. 


0.569 


0.569 


0.099 


0.098 


-0.153 


-0.153 


500. 


0.584 


0.584 


0.117 


0.116 


-0.164 


-0.164 


520. 


0.589 


0.590 


0.138 


0.138 


-0.174 


-0.174 


540. 


0.585 


0.585 


0.163 


0.163 


-0.184 


-0.185 


560. 


0.571 


0.571 


0.192 


0.192 


-0.194 


-0.195 


580. 


0.548 


0.547 


0.226 


0.227 


-0.203 


-0.204 


600. 


0.517 


0.517 


0.267 


0.268 


-0.212 


-0.213 


620. 


0.480 


0.479 


0.315 


0.316 


-0.221 


-0.222 


640. 


0.439 


0.439 


0.370 


0.372 


-0.229 


-0.230 


660. 


0.395 


0.394 


0.433 


0.437 


-0.237 


-0.237 


680. 


0.350 


0.350 


0.496 


0.495 


-0.244 


-0.244 


700. 


0.304 


0.305 


0.540 


0.538 


-0.252 


-0.249 
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TABLE V. Analysis of </>o based on the phase shift data of Ochs and of Estabrooks and Martin 
(E-M), extrapolated to threshold, for fixed values of the scattering length ao, using the parametriza- 
tion of Schenk. 



Data 

J. J CXsKiCXi 


bn 


MJ 




v 2 /d f 

A /U.l. 


a//3 


fn 


v 2 /d f 

A /U.l. 

(imprvd.) 


v 2 /d f 
(stand.) 










ClQ = 


0.20 








Ochs "a" 

V / LAID Oj 


o 393 


-0 0356 

\J .Kit J 'J \ I 


867 8 


n /i 3 










Ochs "b" 


0.348 


-0 0292 


863.3 


10/13 


2.32(4) 


-0 094(6) 


67/63 


690/63 


Ochs "c" 


0.298 


-0.0206 


858.8 


11/13 










E-M "a" 


0.253 


-0.0184 


818.6 


15/15 










E-M "b" 


0.229 


-0.0147 


814.0 


14/15 


1 581(3) 


-0 236(5) 


30/63 


677/63 


E-M "c" 


0.205 


-0.0109 


809.6 


15/15 


















ao = 


0.26 








Ochs "a" 


0.369 


-0.0339 


867.4 


11/13 










Ochs "b" 


0.324 


-0.0274 


863.1 


10/13 


3.50(3) 


-0.154(6) 


48/63 


2328/63 


Ochs "c" 


0.274 


-0.0188 


858.6 


11/13 










E-M "a" 


0.227 


-0.0163 


817.8 


14/15 










E-M "b" 


0.203 


-0.0126 


813.3 


13/15 


2.86(3) 


-0.305(2) 


73/63 


4947/63 


E-M "c" 


0.179 


-0.0087 


808.9 


14/15 
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